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We present a many body analysis of the multi-band Kondo lattice model. The study is then combined with 
the first principles TB-LMTO band structure calculations, in order to investigate the temperature dependent 
correlation effects in the 3p valence bands of the ferromagnetic semiconductor EuS. Some of the physical 
properties of interest like the quasi-particle density of states (Q-DOS), spectral density (SD) and quasi- 
particle band structure (Q-BS) are calculated and discussed. Therewith, we propose a spin resolved ARPES 
of the valence bands of EuS to be performed. 
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1 INTRODUCTION 

The europium chalcogenides EuX, with X=0,S,Se and Te, formed the subject of many studies for more 
than forty years. They are magnetic semiconductors and one of the early motivation was the hope of using 
these materials as magneto-optical memories in the computers and magneto-optical modulators. Moreover, 
the europium monochalcogenides are face centered cubic materials that are ideal examples of Heisenberg 
exchange model and thus formed good exemplary substances for studying magnetism. In such a motivated 
background experiments were performed to study their properties and use them for practical purposes. 
Though apart from various dopings and substitutions, it turned out that the ordering temperature could 
not be raised up to room temperature. Also, these materials were found difficult to prepare. But, despite 
such difficulties, the studies on these materials attracted attention from a more academic point of view 
as it highlighted many interesting physical properties. For example [1, 2], a metal-insulator transition in 
Eu-rich EuO has been observed while red-shift of the absorption edge, Faraday rotation, circular dichroism 
and spin polarization explained magnetic and magneto-optical effects. Such interesting properties have not 
only provided the motivation for theoretical and experimental investigations on these materials, but further 
studies on them permit a better understanding on more complex systems like the Diluted Magnetic Semi- 
conductors (DMS) which are believed to be one of the essential materials behind the technology called 
spintronics [3, 4]. 

The motivation behind such a study is to combine many body theoretical analysis with bandstructure 
calculations for real materials. The goal is to study the temperature dependent electronic correlation ef- 
fects in the valence bands of EuS. A similar study has been carried out earlier for the conduction bands 
of EuS [5] but under a different kind of a theoretical approach. The chalcogenide, EuS, was under a lot 
of experimental and theoretical investigation to study its magnetic and electronic properties, starting from 
early 1960s which reported nuclear magnetic resonance [6] till the latest being on EuS thin films [7] while 
the list remains extensive. 

This paper is organized as follows. In section 2, we formulate the complete multi-band Hamiltonian. 
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A self energy ansatz, reliable for low carrier densities, is used to solve the model Hamiltonian. The tem- 
perature and spin dependent density of states (DOS) are presented for two weakly hybridized bands in a 
simple cubic lattice with coupling strength, hybridization and magnetization being the parameters for the 
model study. The results are then compared with one of the exact limiting cases of the model, namely the 
ferromagnetically saturated semiconductor and are found to be coinciding. In section 3, such a many body 
model is combined with first principles T=0 band-structure calculation, tight binding linear muffin-tin or- 
bital (TB-LMTO), in order to calculate some of the physical properties like quasi-particle density of states 
(Q-DOS), spectral density (SD) and quasi-particle band structure (Q-BS) for realistic ferromagnetic semi- 
conducting material EuS. The final section 4, concludes with the proposition of a spin resolved ARPES 
experiment to be performed in order to examine our results. 



2 THEORY AND MODEL CALCULATION 

The multi-band Kondo lattice model (KLM) Hamiltonian reads as follows; 

H = Hkin + Hi nt (1) 

where 

H kin = T i*/ c LcT C j0<T (2) 

and 

Hint — ~ ^ y ] ( z &Sj c \ aa C iaa + S; t c[ a _ a Ci aa ) (3) 

The term H^n denotes the kinetic energy of the valence band electrons, 3p orbital in case of EuS, with 
c Lct an d c ia° being the fermionic creation and annihilation operators, respectively, at lattice site Ri. The 
latin letters (i,j,...) as subscripts, symbolize the crystal lattice indices while the band indices are depicted 
as superscripts in Greek letters (a, (3,..) and the spin is denoted as cr(= |, J.). Such a notation is used 
throughout the text. The multi-band hopping term, T"? , is connected by Fourier transformation to the free 
Bloch energies e Q/3 (k) 



T f ^E^W e-*'*-*') (4) 

k 



H^t is an intra-atomic exchange interaction term being further split into two subterms. The first de- 
scribes the Ising type interaction between the z-component of the localized and itinerant carrier spins while 
the other comprises spin exchange processes which are responsible for many of the KLM properties. J is 
the exchange coupling strength which we assume to be k-independent and S? refers to the localized spin 
at site Ri 

S? = Sf + iz a S\ ; z T = +1, z t = -1 (5) 

The Hamiltonian in eq. (1) provokes a nontrivial many body problem that cannot be solved exactly. 
Approximations must be considered. We proceed to solve the problem by the equation of motion method 
using double-time retarded Green function [8] 

GZA E ) = ((ci^cL*))e (6) 
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where l,m and are the lattice and band indices respectively. The equation of motion reads as follows 

EGZAE) = hS lm S^ +Y J T??(f} v mrT {E) - J -[TZAE) + F^E)] (7) 

31 

where we use the following notations 

rZA E ) = ((Sid^cl„ a )) E (8a) 

FiZA E ) = ((sr^;cl va )) E (8b) 

We find that eq.(8a) and eq.(8b) are higher order Green functions blocking the direct solution of the 
equation of motion as they can't be decoupled into the forms of the original Green function. But, a rather 
formal solution can be stated as 

G k AE) = hT{(E + iO+)T-e(k)-£ k AE)}- 1 (9) 

where we exclude the band indices by representing the terms in a generalized matrix form on symbolizing 
a hat over it 

G lma {E) = ^E 6 ^) e~ ik < R '- R ^ (10) 

k 

The terms in eq. (9) are explained as follows : / is an identity matrix and e(k) is a hopping matrix with 
the diagonal terms of the matrix exemplifying the intra-band hopping and the off-diagonal terms denoting 
the inter-band hopping. The self energy, S^ ct (E), containing all the influences of the different interactions 
being of fundamental importance can be understood as 

<<[iWw]-;cL CT )> = J2KpAE)G^ m AE) (ii) 

pi 

Now, we are left with a problem of finding a multi-band self energy ansatz, in order to compute the 
Green function matrix and thereby calculate some of the physical quantities of interest like the quasi- 
particle spectral density (SD) 

AUE) = --ImTrG ka (E) (12) 
and the quasi-particle density of states (Q-DOS) 

pAe) = -^Y. a ^ e ) < 13 > 

k 

As observed from our many body theoretical model, we are only interested in the influence of inter- 
band exchange on the valence band states in order to study the electronic correlations and not aimed at 
calculating the magnetic properties via self consistent calculation of the localized magnetization. For this 
purpose, we neglect the exchange interaction between the localized spins, Sf, in our model. Furthermore, 
since we consider fully occupied valence bands, the Coulomb interaction is unimportant and this can further 
reduce the complexity in formulating the self energy ansatz. Actually, in case of realistic calculations of 
EuS, the Coulomb interaction is implicitly taken in the LSDA bandstructure calculations. So, the single 
particle excitation energies contain the Coulomb interactions and need not be taken into account in the self 
energy ansatz. Thus, the criterion required by our multi-band electronic self energy is that it should be 
accurately defined in the low carrier density regime and weak coupling strength, which is the case for the 
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study of magnetic semiconductors. 

It is known that the above statements can be satisfied by the interpolating self energy ansatz (ISA) being 
proposed for a single band model [9]. Such an ansatz is considered to work for all coupling strengths and 
band occupations. Various exactly known limiting cases like the weak and strong coupling behavior, zero 
bandwidth limit [10], etc. are shown to be satisfied by the single-band ISA. Though it is not self-evident 
in case of multi-band situation that these confining cases would hold true, nevertheless we propose to 
replace the single particle Green function for the multi-band case by the respective matrix Green function. 
As seen in theoretical analysis given in Appendix A, the important case of ferromagnetically saturated 
semiconductor gives the same self-energy ansatz as in the single band case (eq.(26) and eq.(27) in Ref. 
9) which strongly supports our ansatz. With the help of such analysis, we put forward for the multi-band 
self-energy (ISA); 

-l 

(14) 
where 

M a = z a {S z )- a ff = S(S + 1) + M a (M a + 1). 
and the bare Green function matrix is defined as : 

d a (E) = 1 y L 

N k (E + iO+)I-e(k) 

It retains the formerly proposed ISA structure in the single band case. As seen in eq. (14), we are 
interested only in the local self energy 

k 

The physical reason for the wave-vector dependence of the self energy is mainly due to the magnon 
energies h(d(k) appearing at finite temperature. But, as stated earlier, we neglect a direct Heisenberg ex- 
change between the localized spins and this can be interpreted as the hui(k) — ► limit. The localized 
magnetization (S z ) shall be considered as an external parameter being responsible for the induced temper- 
ature dependence of the band states. 

In principle, apart from the above discussed proposal one could also consider other alternative ap- 
proaches for the many body analysis like the dynamical mean field theory (DMFT) in order to study the 
correlation effects. But, it is known [11] that the methods of DMFT can not be directly applied to the 
Kondo lattice model. One exception is the case of the classical spin limit, thus removing the quantum 
nature of the spins. While there is another possibility to derive a DMFT for Kondo lattice model, based on 
the fermionization of the localized spin operators as suggested in [12], but it is limited to S=|. Recently, 
in an article [13], the correlation effects were discussed by introducing DMFT based approaches but for a 
single band Kondo lattice model and with an evaluation compared with the exactly solvable limiting cases 
and several other known approximate methods. 

Considering our multi-band electronic self energy ansatz, we proceed to perform a model calculation of 
two weakly hybridized bands on a simple cubic (sc) lattice where the terms in the (2 x 2) hopping matrix, 
e(k) , are taken using the sc Bloch density of states (DOS) in the tight binding approximation [14]. We 
assume the bandwidth (W) of both the bands, i.e intra-band transfer energy to be 1.0 eV while that for 
inter-band to be 0.10 eV. The center of gravity of one of the free Bloch band is chosen to be energy zero 
while for the other it is taken as 0.25 eV. 

In order to study the electronic correlations, one examines the effect of a test electron by creating (or 
annihilating) it in an empty (or filled ) band. In our case, we first annihilate an electron from the completely 
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filled band. Thus, removal of J.(T) electron will create a f(|) hole which will give insight into the study. 
The correlation effects are studied by calculating the electronic self energy. But, apart from electronic 
sub-system, we also have the magnetic sub-system. The exchange coupling between the itinerant electron 
and localized spins adds up to the correlation effects as they produce spin-flip transitions and Ising like 
interactions in addition to the kinetic energy. 

Fig. 1 describes the weak coupling behavior, JS^;W. The calculated density of states (DOS) are shown 
for the parameters : J=0.2 eV, S=1.5 and different values of magnetization, 1 ^r 1 =1.0, 2/3, 1/3 and 0.0 
which are represented as thick solid line, dotted line, broken line and double-dotted broken line, respec- 
tively. The exact solution, ferromagnetically saturated semiconductor, is represented as stars and it is found 
to coincide with results of (5 Z )=S. 
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Fig. 1 The spin up (<r=|)and down (<7=J.) DOS, as func- 
tion of energy, are calculated with J=0.2 eV, S=1.5 for var- 

(s z ) 

ious values of magnetization ^-^-=1.0, 2/3, 1/3 and 0.0 
represented as thick solid line, dotted line, broken line and 
double-dotted broken line, respectively. The case of fer- 
romagnetically saturated semiconductor is represented as 

(s z ) 

stars coinciding with -^-=1.0. 



J=1.0 eV, S=1.5 



O 

a 




Fig. 2 The same as Fig. 1 except for J=0.5. 



Fig. 3 The same as Fig. 1 except for J=1.0. 



At T=0 K and fully occupied bands, the removal of a j electron (i.e. creation of a f hole) produces 
a stable quasi-particle since the | hole has no chance to flip its spin with the ferromagnetically saturated 
spin (ff) sub-system. The imaginary part of the self energy vanishes indicating infinite lifetimes. The sit- 
uation is different for the removal of an | electron (i.e, creation of J, hole) as the J, hole can exchange its 
spin with the ferromagnetically saturated spin sub-system and have finite lifetime. It can emit a magnon 
and become a f hole provided there exist f hole states which can be occupied after spin-flip process. These 
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are the scattering states which occupy the same energy regions as of the f spectrum. The scattering states 
become more distinct with increasing coupling strength as seen in Fig. 2 and Fig. 3. The other possibility 
for the J, hole to exchange its spin could be by repeated emission and absorption of magnons, thus forming 
a quasi-particle called magnetic polaron. 

However, at finite temperatures, the spin sub-system is no longer perfectly aligned. There are magnons 
in the system that can be absorbed by the itinerant charge carriers. As seen from the figures, at finite 
temperatures the spectral weight gets redistributed due to spin flip term in the exchange interaction with 
deformations in the density of states. In the limit T — > T c ((S z ) — > 0) the induced spin asymmetry is 
removed. 



3 PHYSICAL PROPERTIES : EuS 



In this section, we study the electronic correlations in the 3p valence bands of ferromagnetic semiconductor 
EuS by combining the multi-band self energy ansatz for the many body part along with the first principles 
TB-LMTO bandstructure calculations. 

EuS crystallizes in a rocksalt structure with lattice constant, a=5.95 A. Each Eu 2+ ion has twelve 
nearest and six next nearest Eu-neighbors and they occupy lattice sites of a fee structure. The magnetism 
in EuS is due to the half-filled 4/ shell of Eu 2+ and its magnetic properties are well described using the 
Heisenberg model with the ferromagnetic transition temperature, T c =16.54 K. In this material, the nearest 
neighbor exchange interaction is ferromagnetic while the next nearest neighbor is anti-ferromagnetic. 

In order to have the single particle excitation energies, i.e the hopping matrices, and all the interactions 
which are not directly included in our Hamiltonian, we perform the TB-LMTO bandstructure calculations 
using the program of Anderson [15, 16]. 

Fig. 4 shows the spin-dependent bandstructure (LSDA) of the valence bands of EuS. In TB-LMTO, the 




Fig. 4 LDA bandstructure of 3p valence bands of EuS. Full lines are spin 
| and dotted ones are spin J.. 

original Hamlitonian is transformed to a tight-binding Hamlitonian containing only the nearest neighbor 
correlations. The evaluation is restricted only to 3p bands. There are some difficulties which arise due to 
the strongly localized character of 4/ levels and in order to avoid such difficulties, we consider the seven 4/ 
electrons as core electrons, since these electrons enter our model only as localized moments. 

The LDA-density of states is displayed in Fig. 5. Assuming that the LDA treatment of the ferromag- 
netism is quite compatible with the S toner (mean field) picture, the T=0 splitting amounts to AE=JS. On 
taking the centers of gravity of the DOS for both the spins and along with the above assumption, the ex- 
change splitting amounts to AE=0.1512 eV which results in the value for the exchange coupling strength 
as J=0.0432 eV where the spin value, S=3.5. One can also calculate the exchange coupling by taking the 
energy difference of the upper edge shift between both the spins which gives J=0.0253 eV or calculate the 
mean value of the spin polarized energy difference taken along each k-points in the first Brillouin zone 
which gives J=0.0529 eV. It was found that the coupling strength obtained by center of gravity approach, 
reflects more correlation effects as was shown in the previous work [5] for the conduction bands of EuS. 

The single particle output obtained from the bandstructure calculations are in the form of Hamiltonian 
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Fig. 5 Spin-dependent density of states of 3p bands of EuS. Using 
the center of gravity of the bands, the exchange splitting amounts 
to AE=0.1512eV. 




and overlap matrix, posing a generalized eigenvalue problem to be solved. In order to employ such ma- 
trices as an input for the many body calculations, one has to perform a Cholesky decomposition so as to 
reduce the generalized problem to an eigenvalue problem. Using such an approach, one can factorize the 
Hamiltonian matrices and then obtain new Hamiltonian matrices using overlap matrices, which can be used 
directly as hopping integrals, e(k). We use the multi-band self energy ansatz eq. (14) to compute the Green 
function eq. (6) and therewith calculate the spectral densities eq. (12) and densities of states eq. (13). 
Fig. 6 represents the quasi-particle band-structure for some high-symmetry directions in the first Bril- 
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Fig. 6 Spin-dependent quasi-particle bandstructure of 3p valence bands of EuS for different values of magnetization 
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louin zone. The degree of blackening measures the magnitude of the spectral function. Correlation effects 
are clearly visible along many parts, especially along the center of the valence band. As seen in the f 
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spectrum even at T=0 K (■^j =1), parts of the dispersions are washed out showing lifetime effects due 
to correlation in terms of magnon emission and re-absorption with simultaneous spin-flips as explained 

earlier. But at the same temperature, the j spectrum shows sharp lines indicating infinite lifetime as ex- 

(s z ) 

pected. While at finite temperature (-^^=0.5), the finite lifetimes appear in case of both the spin spectra 

as the ferromagnetically saturated spin system is not perfectly aligned giving rise to magnons which can 

(s z ) 

be absorbed by the itinerant charge carriers. At transition temperature ( 1 -j^=0.0), the spin asymmetry is 
removed. Strong temperature dependent correlation effects are seen mainly along the W-L directions. 

In order to have a closer look at the correlations, the spin and k-dependent spectral densities are plotted 
at some high-symmetry points (W,L,X) in the first Brillouin zone. It can be compared with the angle- 
resolved photoemission experiments and for this purpose we propose that an experiment be performed. 
As seen in Fig. 7, well defined quasi-particle peaks appear with an additional spin split below T c . The 
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Fig. 7 Spectral density at X - point. 



exchange splitting collapses for T — > T c and the quasi-particle damping increases with increasing tempera- 
ture. Similar explanations hold for the other high-symmetry points (Fig. 8 and Fig. 9). The correlations are 
clearly observed at the W-point. With increasing temperature, there's a strong damping as seen along the 
energy spectrum. Interestingly, the same W-point remains the point of discussion in case of the conduction 
band calculations of EuS T51. 



mi 




<s >/s = 
<s 7 ">/s = 
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Fig. 8 Spectral density at W - point. 



Fig. 9 Spectral density at L - point. 



© 2003 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



pss header will be provided by the publisher 



9 



The quasi-particle density of states (Q-DOS) is plotted in Fig. 10. As observed, due to weaker ex- 
change coupling, the appearance of polaronlike quasi-particle branches are less likely in the real material 
calculation. The temperature influence on the spectrum is seen in terms of strong deformations, mostly 
at the center, and in terms of shifts. While the k-dependent spectral densities refers to the angle resolved 
photoemission, the Q-DOS can be considered as the angle integrated part. 



<S*>/S=1.00 

<S*>/S = 0.75 



- <S Z >/S = 0.50 
- - <S Z >/S = 0.25 




-6 -5 -4 

Energy 



Fig. 10 Quasi-particle density of states (Q-DOS) of the 3p valence bands of EuS. 



4 SUMMARY and CONCLUSION 

In this paper, we have studied the effects of electronic correlation on the 3p valence bands of the ferromag- 
netic semiconductor EuS. We started with multi-band Kondo lattice Hamiltonian which is sufficient enough 
to explain the necessary physics. Using Green function approach, we try to evaluate the Hamiltonian. We 
put forward an Interpolating Self energy Ansatz (ISA) for the multi-band situation and obtain the solution, 
i.e the required Green function, to predict some physical properties. In order to test the theory, we perform 
model calculations and interesting features are highlighted. For two weakly hybridized bands assuming 
free Bloch densities of states on a simple cubic lattice, polaron-like and scattering states are observed for 
intermediate and strong coupling strengths. Strong temperature dependent correlation effects are noticed 
in the form of deformations and shifts in the energy spectrum. A non-trivial exactly solvable limiting case 
of the model is correctly reproduced thus confirming the proposed multi-band self energy ansatz. 

Such a self energy ansatz, which amounts to the many body interactions, is later on combined with the 
first principles TB-LMTO bandstructure calculations for calculating the properties of real materials. The 
bandstructure output in the form of hamiltonian and overlap matrices of the 3p valence band electrons of 
EuS, serve as an input in the form of hopping integrals for the many body calculations. The exchange cou- 
pling strength is obtained from the LDA bandstructure calculations as a result of energy difference in the 
center of gravities of the spin polarized density of states. With these parameter inputs, physical properties 
like quasi-particle spectral densities (Q-SD) and densities of states (Q-DOS) are estimated for different 
values of magnetization. One can also obtain other interesting physical quantities like the lifetime or the 
effective mass of the quasi-particles but it is not planned for the present article. While, in one of the pre- 
vious work [17], emphasis was laid upon calculating the spin flip probabilities which exhibit the lifetime 
effects. These quantities were calculated using the time-dependent spectral density which were used to 
express the spin polarization of the field emitted electrons thus highlighting the issue of spin filter effects 
of the europium chalcogenides. Eventually, as a consequence of our present calculations, the temperature 
dependent correlation effects are observed at the center of band and mainly in the W-L direction. For the 
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sake of comparison, we propose and spin-dependent ARPES to be performed on the valence bands of EuS. 



Acknowledgements One of us (A.S) would like to thank Dr.W.Miiller for his introduction to the TB-LMTO program 
and for helpful discussions. 

A Case of ferromagnetically saturated semiconductor 

In the main text of this paper, we mentioned one limiting approach of the multi-band model namely the 
ferromagnetically saturated semiconductor which exhibits the quasi-particle called magnetic polaron. It is 
that quasi-particle where a bare electron (or hole) is dressed by a virtual cloud of magnons. In this section, 
we present the analytical result of it. We are interested in finding a solution for the model at T=0 with 
ferromagnetically saturated spin sub-system ((S z ) = S) and for fully occupied bands and a=|,|. 
In such situation, the Ising-like higher order Green function eq. (8a) becomes; 

rZAE) = SGZAE) (15) 

Now, let us consider : a) u = f and b) a = J. and try to evaluate the self energy in both the cases, 
a) For a = f, the equation of motion for the spin-flip function 

FZAE) = ((Src hl - c ,;cl va )) E 
reads as follows, 

EFZ l (E) = {{[Src l , l ,HUcl^)) E 
which takes the form as given below yielding higher order Green functions, 

pi 

+£{(£«*r i c Lv]))e) - 2((Sf4 Ai Q AT c^ i ;c^ T ))_ E }] 

A <T 



The higher order Green functions resulting from the spin-flip equation of motion can be evaluated 
exactly for the case of fully occupied bands and at T=0, 

({SfS+c w ;cl^)) E =0 (16a) 
{{STSfci^cl^E = SF^(E) (16b) 

((^r( c L T C(A T - c\\i c im)ci»1\ cLt» B = 5 ^ F ^{E) (16c) 
{{S!c\ Xi c mCllil ;cl^)) E = -S^SGfyE) (16d) 
On substituting the eq. (16) in the spin-flip equation of motion and taking the Fourier transform we get; 

JO 7 7 

(E~ -f)F k ^(E) =^^(k.q)F k ^ qt (£) - JL £ (E) ^=G%{E) 

and in the matrix form, it can be written as follows 

7C 7 7C 

[(E - ^k-q)]F k , k . q , qT (£;) = - — F*mME) - ^ kT (£) 

© 2003 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



pss header will be provided by the publisher 



11 



Let us now define the effective Green function, 

JS_ 1_ 

2 N 



q 

which upon substituting in the matrix form of the spin-flip equation of motion implies, 

1 ^ ^ TO ^ 7^ 70 ^ 

^=^F Kk . q ^(E) = -JSGf{E-—)[I+-Gf{E-—)]-^{E) (17) 

Now , if we consider the Fourier transform of eq. (7) and substitute eq. (17) and eq. (15) in it , we get 
the result for the Green function as 

[ { E+^)T-e{k)]G k ,(E) = hl+ ^Gf(E-^)[T+ J -Gf {E - J -l)]-^{E) 

i.e, 

G kT (£) = W[(£ + i0 +^^ (18) 
On comparing eq. (9) and eq. (18) we get the self energy ansatz as 

^(E) = T+ '^Gf{E - ^)[/+ ldf{E - (19) 
b) For cr = l, the spin-flip function vanishes; 

The equation of motion for the Green function, i.e, eq. (7) reduces to, 

pi 

And upon taking Fourier transform we have, 

TO 

(E —)G^(E) = U, v + J2 ^(*)G J h(E) 

7 

which can be written in matrix form 

G kl (E) = hI[(E - — + iO+)I - ^k)]- 1 (20) 

Thus, the self energy becomes; 

^ tq ^ 

Zi(E) = —I (21) 

As seen, eq. (19) and eq. (21) give the spin-dependent form of the self energy which for the case of fully 
occupied bands and ferromagnetically saturated spin sub-system can be written down as follows; 

^(E) = -^pLT+ j2S ^ + 1) Glff(E ^)[T+ J -&Jf{E - (22) 
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